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ABSTRACT 

New observations of SN 1980K made with the VLA at 20 and 6 cm from 
1994 April through 1996 October show that the supernova (SN) has undergone 
a significant change in its radio emission evolution, dropping by a factor of ~ 2 
below the flux density S oc t -0 73 power-law decline with time t observed earlier. 
However, although S at all observed frequencies has decreased significantly, its 
current spectral index of a — —0.42 ± 0.15 (S oc u +a ) is consistent with the 
previous spectral index of a — — 0.60t[j;o7. 

It is suggested that this decrease in emission may be due to the SN shock 
entering a new region of the circumstellar material which has a lower density 
than that expected for a constant speed (w), constant mass-loss rate (M) wind 
from the progenitor. If such an interpretation is correct, the difference in wind 
and shock speeds appears to indicate a significant evolution in the mass-loss 
history of the SN progenitor ~ 10 4 years before explosion, with a change in 
circumstellar density (oc M/w) occurring over a time span of ^ 4 kyr. Such 
features could be explained in terms of a fast "blue-loop" evolutionary phase 
of a relatively massive pre-SN progenitor star. If so, we may, for the first time, 
provide a stringent constraint on the mass of the SN progenitor based solely on 
the SN's radio emission. 

Subject headings: radio continuum: stars, supernovae: individual (SN 1980K) 
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1. Introduction 

SN 1980K [RA(1950.0)= 20 h 34 m 26!68 ± OfOl; Dec(1950.0)= +59°55'56'. / 5 ± 0"2] in NGC 
6946 was discovered in 1980 October by Wild (1980), and an initial, unsuccessful attempt 
was made on 1980 November 3 to detect it at 6 cm with the Very Large Array (VLA).Q 
However, only 35 days after optical maximum SN 1980K was detected at 6 cm on 1980 
December 4 and has been regularly monitored at 20, 6, and occasionally 2 cm since then. 
Results from 1980 November through 1984 August are presented in Weiler et al. (1986), 
and results from 1984 November through 1990 December are presented in Weiler et al. 
(1992). By 1990 December the flux density had fallen to S < 0.3 mJy at 20 and 6 cm, which 
made it difficult to observe with the VLA using typical, short "snapshot" observations, and 
regular monitoring was terminated. 

However, optical imaging and spectra taken after 1990 indicated that SN 1980K was 
still a detectable, albeit relatively constant optical luminosity source, particularly at Ha 
(Leibundgut 1994). Late-time optical and radio emission are predicted to be correlated 
(Chevalier & Fransson 1994), and in several examples are indeed observed to be related (see, 
e.g., Montes et al. 1997, Van Dyk et al. 1998 and references therein). Therefore, deeper and 
more sensitive VLA observations were made in 1994, 1995, and 1996 at both 20 and 6 cm, 
to attempt to obtain new measurements of SN 1980K below our previous sensitivity limits. 
Interestingly, these new data indicate that SN 1980K has, in fact, not continued along 
its previously observed power-law decline, but has dropped sharply in radio flux density 
in the interval between 1990 and 1994. Additionally, recent measurements (Fesen 1998) 
indicate that SN 1980K has also faded in the optical over the last several years, possibly 

4 The VLA is operated by the National Radio Astronomy Observatory (NRAO) of 
Associated Universities, Inc., under a cooperative agreement with the National Science 
Foundation. 
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starting with a 20% decline by 1994 as reported by Fesen, Hurford, & Matonick (1995). 
The new radio measurements are presented here, and they indicate a significant change in 
the evolution of the radio emission of SN 1980K from earlier, model-based expectations. 

2. New Observations 

The new observations of SN 1980K were made with the VLA at 20 cm (1.46 GHz) and 6 
cm (4.89 GHz) on 1994 April 10, 1995 June 22 (20 cm only), and 1996 October 13-14, when 
the VLA was in the "A" configuration, to reduce the confusion with background emission 
from the parent galaxy NGC 6946. Each observation was accompanied by measurement 
at each frequency of a "secondary," possibly variable, calibrator (2021+614), with a 
well-known position, making it suitable as a phase reference. To establish an absolute flux 
density scale, a "primary" calibrator (3C 286, except on 1995 June 22, when 3C 48 was 
used) was also observed and used to calibrate the amplitude of 2021+614 for each observing 
session. Additional observations of SN 1980K on 1994 April 01 and 1994 June 16 were 
kindly provided by C. Lacey (1997; see also Lacey, Duric, & Goss 1997). 

2.1. Calibration 

The VLA is described in a number of publications (e.g., Thompson et al. 1980; 
Hjellming & Bignell 1982; Napier, Thompson, & Ekers 1983); the general procedures for 
RSN observations, with the sources of possible error, are discussed in Weiler et al. (1986). 
The primary calibrators employed here, 3C 48 for the 1995 June observation and 3C 286 
for the rest, were assumed to be constant in flux density with time and to have the flux 
densities at 20 cm and 6 cm of S 20 cm (3C 48) = 15.97 Jy, cm (3C 286) = 14.75 Jy, and 
S 6 cm (3C 286) = 7.49 Jy. 
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The secondary calibrator, 2021+614, was used as the phase (position) 
[RA(1950.0)=20 h 21 m 13 f 5300, Dec(1950.0)=61°27'18"156] and amplitude (flux density) 
reference for the observations of SN 1980K. Its flux density (measured against 3C 48 or 3C 
286) at 20 cm and 6 cm is presented for each observing date in Table fjj and its flux density 
evolution with time since 1980 is shown in Figure |l|. The data were processed in a standard 
manner using the Astronomical Image Processing System (AIPS) provided by NRAO. 

2.2. Error Estimates 

Error estimates are derived as in Weiler et al. (1991). That is, the rms "map" error, 
o"o, in an image only provides a lower limit to the total error, so that we add quadratically 
an error of 5% of the measured flux density to account for inaccuracies in the flux density 
calibration and any absolute scale errors in the primary calibrators. The the final errors, 
a - /, are taken as 

a) = (0.05S ) 2 + cx 2 , (1) 

where So is the observed flux density for SN 1980K and o"o is the observed rms map error 
measured outside of any obvious regions of emission. 

3. Results 

Since the previous publication of SN 1980K data in Weiler et al. (1992), we have added 
three new VLA measurements at 20 cm, and two at 6 cm. 

NGC 6946 was also surveyed with the VLA by Lacey et al. (1997) for compact sources 
on 1994 April 1 at 20 cm in the "A" configuration, and on 1994 June 16 at 6 cm in the "B" 
configuration. Although outside of their 9' x 9' region of interest centered on the nucleus of 
NGC 6946, SN 1980K was detected in their survey at 20 cm. Since the SN was ~ 5.5' from 
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their field center and suffered primary beam attenuation, their measurement of flux density 
was corrected by 6.9% with the standard AIPS task PBCOR. Their 6 cm observation 
suffered even more severe attenuation at the location of SN 1980K, and only an upper limit 
on the 6 cm flux density could be established for the SN. 

Twelve revised measurements^ between 1980 November 3 and 1990 December 9 and 
the six new measurements plus one upper limit between 1994 April 1 and 1996 October 14 
are presented in Table |[ In Figure |] we plot the entire time evolution of the flux density 
for SN 1980K at both 20 and 6 cm, from 1980 November 3 to 1996 October 14. This 
includes the new (and revised) results, together with earlier measurements from Weiler et 
al. (1986) and Weiler et al. (1992). The 6 cm upper limit from Lacey (1997) for 1994 June 
16 is consistent with our measurement on 1994 April 10. The 20 cm measurement by Lacey 
for 1994 April 1 is about 1.5a greater than our 20 cm measurement of 1994 April 10, but 
the reported uncertainties overlap. 

4. Parameterized Model 

Weiler et al. (1986) discuss the common properties of radio-emitting supernovae 
(SNe), known as radio SNe (RSNe), including nonthermal synchrotron emission with 
high brightness temperature, turn-on delay at longer wavelengths, power-law decline after 

5 We include in Table ^] eight previously published measurements which were originally 
listed in Weiler et al. (1986, 1992) as < 3a detections. We now feel a conservative approach to 
consider these data as 3a upper limits is more appropriate. Additionally, we have remeasured 
four previously-published results and list these revised values in Table ||. As can be seen 
from a comparison of Table § in this paper with Table 4 in Weiler et al. (1992), the effect of 
these changes on the best-fit model parameters is relatively small. 
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maximum with index /3, and spectral index a asymptotically decreasing to an optically 
thin value. Weiler et al. (1986) have also shown that the "mini-shell" model of Chevalier 
(1982a,b) adequately describes known Type II RSNe. In this model, the relativistic 
electrons and enhanced magnetic fields necessary for synchrotron emission are generated by 
the SN shock interacting with a relatively high-density ionized circumstellar envelope. This 
dense cocoon is presumed to have been established by a high mass-loss rate (M > 
M Q yr -1 ), low speed (u> w ind ~ 10 km s _1 ) wind from a red supergiant (RSG) SN 
progenitor which was ionized and heated by the initial SN UV/X-ray flash. The rapid rise in 
radio flux density results from the shock overtaking progressively more of the wind matter, 
leaving less of it along the line-of-sight to absorb the emission from the shock region. 

Following Weiler et al. (1986), we adopt the parameterized model: 



and with K\ and K 2 corresponding, formally, to the flux density in mJy and the uniform 
external absorption, respectively, at 5 GHz one day after the explosion, £ - The term e~ T 
describes the attenuation of a local medium that uniformly covers the emitting source 
("uniform external absorption"). This absorbing medium is assumed to be purely thermal, 
ionized hydrogen with frequency dependence u~ 21 . The parameter 5 describes the time 
dependence of the optical depth for this local, uniform medium. For an undecelerated SN 
shock, 5 = —3 is appropriate (Chevalier 1982a). 

This parameterization has been found generally applicable to Type II SNe, such as 
SN 1979C (Weiler et al. 1986, 1991) and SN 1980K (Weiler et al. 1986, 1992), with values 
of 5 close to the undecelerated value (5 = -3.12 and -2.68 for SN 1979C and SN 1980K, 
respectively; see Table |3|). 





where 




(3) 
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The best-fit model parameters are listed in Table [| only using data through 1989 July 
23. This fit is slightly different from that of Weiler et al. (1992), due to improved fitting 
software and reclassification and remeasurement of some data points (see Footnote 0), but 
the difference is not deemed to be significant. 

The errors in the fitting parameters in Table |^ were estimated using a "bootstrap" 
procedure (Press et al. 1992). Bootstrap procedures use the actual data sets to generate 
thousands of synthetic data sets that have the same number of data points, but some 
fraction of the data is replaced by duplicated original points. The fitting parameters are 
then estimated for these synthetic data sets using the same algorithms that are used to 
determine the parameters from the actual data. The ensemble of parameter fits is then used 
to estimate errors for the parameters by examining number distributions for the parameter 
in question. The errors we present correspond to the 15.85% and 84.15% points in the 
distribution of the occurrence of parameters from the synthetic data. 

Figure ^| shows all of the data available through 1996 October 14, with the best-fit 
model up to day 3200 shown for 6 cm (solid line) and 20 cm (dashed line). After day 4900 
we extrapolate this model with the same decline rate /3 and same spectral index a, but with 
a discontinuity of 55% introduced between day 3200 and day 4900. Note that the single 
detection on day 3692 (1990 December 9), and several low upper limits at 6 cm, imply that 
the drop-off actually started sometime between the last 6 cm detection on day 3188 (1989 
July 23) and day 3692. Also note that the extrapolated model light curves after day 4900, 
while showing the preservation of the spectral index (see Figure |3[) and consistent with 
the available data, are not considered to have any predictive value, because of the large 
uncertainty in the data. The data are also consistent with a constant flux density (i.e., with 
/3 = 0) after day 4900. 

Figure |3| shows the evolution of the spectral index, af, between 6 and 20 cm; the solid 
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line is derived from the model fit in Table 0. 

5. Discussion 

Examination of Figure Q shows that radio observations of SN 1980K, at least through 
day 3200, are consistent with the tenets of the Chevalier (1982a,b) model, as described by 
Weiler et al. (1992). However, sometime between day ~ 3200 and day ~ 4900 the radio 
emission of SN 1980K changed dramatically, dropping by a factor of ~ 2 below the values 
expected from the previously best-fit model. This roll-off is even roughly described by the 
single 6 cm detection on day 3692 (1990 December 9) and the several 6 cm upper limits 
around that time. Interestingly, however, the large change in flux density has apparently 
not affected the spectral index evolution shown in Figure || The best-fit post-decrease 
spectral index, «g°(day > 4900) = —0.42 ± 0.15, is consistent to within the errors with 
the pre- decrease spectral index, c^°(day < 3200) = -0.60lg$. While such a sharp drop 
could imply a change in emission mechanism, emission efficiency, or shock transition to 
a non-radiative phase, we interpret the apparently constant spectral index to imply that 
the basic shock-driven emission mechanism has remained unchanged, and that the flux 
density drop is most likely due to a deviation of the average circumstellar medium (CSM) 
density behavior from the standard r~ 2 law, expected for a pre-SN RSG wind with constant 
mass- loss rate (M) and constant speed (w). 

/,y\ (7-7+12m)/4 

Since the radio luminosity of a SN is L ra dio oc J (7 = —2a + 1 is the 

relativistic electron energy index and m = —5/3, where the shock radius evolves as R oc t m , 
and 5 is defined in §4; Chevalier 1982a,b), it is related to the average circumstellar density 
(Pcsm oc M/w), a measure of the deviation of the radio emission from the standard model 
may provide, assuming all other aspects of the emission mechanism remain unchanged, an 
indication of a deviation of the circumstellar density from the canonical r~ 2 law appropriate 
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for a constant mass- loss rate (M), constant speed (w) pre-SN stellar wind. Thus, for 
SN 1980K with L radio oc (M/u>) L53 , a flux density decrease by a factor of ~ 2 implies, 
for a constant emission mechanism, a decrease in pcsm by a factor of ~ 1.6. [N.B.: The 
expressions we use from the Chevalier model assume a power-law CSM density distribution. 
They predict that a very small CS density leads to a very small radio luminosity. However, 
if there is a sudden drop in the density, the previously accelerated electrons will continue 
to radiate, and the radio luminosity will only decrease due to adiabatic expansion. This 
implies that the density changes which we obtain are, strictly speaking, lower limits, since 
some of the radio luminosity must be coming from electrons that were accelerated before 
the shock entered the less dense medium.] 

Weiler et al. (1986) predicted the possibility of a sharp decrease in the flux density for 
RSNe when the shock reaches the edge of the high-density, low-speed stellar wind CSM 
established during the RSG phase of the pre-SN star, and such a sharp decrease in flux 
density has also been observed in SN 1957D at 6 and 20 cm (Cowan, Roberts, & Branch 
1994), accompanied by a steep decline at optical wavelengths (Long, Winkler, & Blair 
1992). Unfortunately, since the SN 1957D radio measurements span only two epochs with 
the observations occurring ~26 and ~33 years after discovery, we have no knowledge of 
the prior radio behavior (spectral index evolution and temporal history) of SN 1957D, and 
such a sharp drop in flux density has not been documented in any other supernova with 
well-observed radio light curves. In at least the two cases of SN 1987A (Turtle et al. 1990; 
Staveley-Smith et al. 1992, 1993, 1995; Ball et al. 1995; Gaensler et al. 1997) and SN 1979C 
(Montes et al. 1998), changes in the radio flux density evolution at ages of a few years imply 
that the shock is encountering a denser, rather than less dense, CSM. 

Our suggested change in the progenitor's mass-loss rate (or wind speed, or both) 
would have occurred at a time before explosion of about t c han g c ~ 11 ± 2(f shock /w) years 
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or, assuming f s hock ~ 10 4 km s -1 and u> w i n d ~ 10 km s _1 , £ c hangc ~ 11,000 years before 
the SN. (Note that m = —5/3 ~ 0.9 implies a small deceleration of the shock with time, 
since the shock radius R oc t m , but given other uncertainties in this estimate, we will 
not correct for that here.) The change between implied mass-loss rates (or wind speeds, 
or both) occurred within At change <^ 4(t> shock /w) years, or At change <^ 4000 years with the 
assumptions made above. (On day 3200 SN 1980K still appeared to be declining normally; 
by day 3700 it was already suspiciously below the model prediction, and by day 4900 it 
was observed to have declined significantly below model predictions.) The minimum time 
scale of a high-to-low density transition is set by light travel time and corresponds to 
At min ~ 2R/c = 2v s hock/ct ~ 0.09t, so that the drop over ~ 4 years occurring ~ 11 years 
after explosion does not violate this At min ^ 1 year limit. 

Both the pre-1990 evolutionary phase and the 1990-to-1994 change-of-phase are short, 
compared to the lifetime of typical RSG progenitor stars (5-1 Myr for stellar masses 
7-15 Mq, respectively). They are also shorter than the evolutionary time since the last 
"blue-loop" He-burning episode (several xlO 5 years) for stars with masses of ~ 7-12 Mq 
for solar metallicity (and higher masses for lower metallicities; see Brocato & Castellani 
1993; Langer & Maeder 1995). The fast drop observed for the SN 1980K radio emission, 
if due to a density transition, requires a highly-structured CSM from a rapid pre-SN 
stellar evolution and variations of M/w by factors of at least 1.5-2 over time intervals 
as short as 10,000 years or less. Such rapid phase changes have not been explored fully by 
stellar evolution models, but may be similar to the well-known blue-loops occurring in the 
He-burning evolutionary phase, with much shorter time scales. It is conceivable that these 
short excursions are the final phases of evolution for massive stars, as they pass through 
successive nuclear burnings (C-burning, O-burning, etc.), which are of increasingly shorter 
duration and may cause the star to quickly move in the Hertzsprung-Russell diagram, from 
red to blue (or less red), and back again. In this context, it is interesting to note that for 
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SN 1987A the final blue supergiant phase, following a RSG phase, is estimated to have 
lasted < 20, 000 years. It thus appears that the final stages of evolution of massive stars 
may be more complex than usually understood. 

The previously-noted change in the radio evolution of SN 1987A and SN 1979C, 
together with this change in SN 1980K, now provides three examples of SNe where possible 
changes in the CSM density structure have been observed and imply significant changes in 
stellar evolution on time scales of only a few thousand years for the progenitors of these 
SNe. However, since the observed radio characteristics (and, therefore, presumed CSM 
structure) are different in each case, there is apparently a range of evolutionary tracks for 
the massive stars which are SN progenitors. 



6. Conclusions 

At an age of ~ 10 years, SN 1980K has apparently, and unexpectedly, entered a new 
phase of radio evolution, characterized by a large decrease in the flux density at both 6 
and 20 cm below that expected from the model describing its previous power-law decline. 
The spectral index, however, is apparently unchanged and is consistent with its earlier 
value, indicating that the emission mechanism is likely the same. We interpret this radio 
decrease to imply that the circumstellar density has abruptly changed and the SN shock is 
now entering a CSM region ^1.6 times less dense than the RSG wind-established CSM 
previously encountered by the SN shock. This is consistent with models for an ~ 9Mq 
star that undergoes a "blue-loop" evolutionary episode before entering its final RSG phase. 
It is interesting to note that changes in the radio emission, and, by implication, in the 
circumstellar density structure, have also been observed on time scales of a few years (a 
few xlO 3 years, in the time frame of the SN progenitor) for SN 1987A (Turtle et al. 1990; 
Staveley-Smith et al. 1992, 1993, 1995; Ball et al. 1995; Gaensler et al. 1997) and SN 1979C 
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(Montes et al. 1998). 

Although the radio emission from SN 1980K is now very faint 200/iJy), we intend 
to continue monitoring the SN with the VLA, to attempt to determine the reality, extent, 
and the radial profile for this proposed new CSM region. This may allow us to establish 
the properties and duration of this phase of pre-SN stellar evolution, and may, for the first 
time, provide a stringent constraint on the mass of the SN progenitor based solely on the 
SN's radio emission. 
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Days since 1980 October 30 



Fig. 1. — Flux density evolution of the secondary calibrator, 2021+614, with 6 cm (filled 
triangles) and 20 cm (filled diamonds) data shown. Smoothed curves (solid for 6 cm and 
dashed for 20 cm) have been added to guide the eye; the curves were not used in the 
calibration procedure. 
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Fig. 2.— Radio light curves for SN 1980K in NGC 6946. The data at wavelengths 20 
cm (crosses) and 6 cm (open circles) are shown together. The data represent 16 years of 
observation from 1980 November 3 through 1996 October 14, including new and revised 
observations presented in this paper and previous measurements from Weiler et al. (1986) 
and Weiler et al. (1992), as well as measurements by Lacey (1997). The solid and dashed 
curves up to day 3200 represent the best-fit model light curves at 6 and 20 cm, respectively. 
The curves after day 4900 are the extrapolation of the model appropriate up to day 3200 
with the same decline rate f3 and same spectral index a, but with a discontinuity of 55% 
introduced between day 3200 and day 4900. 
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Fig. 3. — The radio spectral index, «g , between 6 and 20 cm for SN 1980K. Notice that, 
even after the the sharp flux density decrease after day ~ 3200 (see Fig. ^), the spectral 
index has apparently remained constant and is consistent, to within the errors, with the 
spectral index from earlier epochs. The solid curve is the spectral index as calculated from 
the model parameters in Table |3], using only data from before (< day 3200) the sharp flux 
density decline seen in Fig. |2|. 
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Table 1. Measured Flux Density Values for the Secondary Calibrator 2021+614 



Observation Primary Flux Density 
Date Calibrator S20 Sq 

(Jy) (Jy) 



1994 Apr 10 3C 286 2.122 2.936 

1995 Jun 22 3C 48 2.124 

1996 Oct 13 3C 286 2.089 

1996 Oct 14 3C 286 • • • 3.086 
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Table 2. New and Revised Flux Density Measurements for SN 1980K a 
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1980 Nov 03 c 


4 


A 


4.885 


<0.300 


0.100 


1981 Feb 04 c 


97 


A 


1.465 


<0.360 


0.120 


lyoz iNOV 11 


( DO 




1.400 


<0.210 


0.070 


1982 Dec 30 c 


791 


D 


1.465 


<0.210 


0.070 


1983 Jan 20 c 


812 


D 


1.465 


<0.210 


0.070 


1983 Feb 18 c 


841 


C 


1.465 


<0.210 


0.070 


1986 Dec 23 c 


2245 


C 


4.885 


<0.450 


0.150 


1987 Mar 16 c 


2328 


D 


4.885 


<0.360 


0.120 


1990 Feb 15 d ' c 


3395 


A 


4.885 


<0.222 


0.074 


1990 May 29 d ' c 


3498 


A 


4.885 


<0.180 


0.060 


1990 Aug 27 d - c 


3588 


B 


4.885 


<0.195 


0.065 


1990 Dec 09 c 


3692 


C 


4.885 


0.181 


0.043 
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Table 2 — Continued 



Observation 


Time Since 


VLA 


Frequency 


Flux Density 


Date 


Optical Maximum 13 


Configuration 


V 


S 


a 




(days) 




(GHz) 


(mJy) 


(mJy) 


1994 Apr 01 f 


4901 


A 


1.450 


0.200 


0.041 


1994 Apr 10 


4910 


A 


1.425 


0.139 


0.041 


1994 Apr 10 


4910 


A 


4.860 


0.104 


0.033 


1994 Jun 16 d ' f 


4977 


B 


4.860 


<0.210 


0.070 


1995 Jun 22 


5348 


A 


1.425 


0.190 


0.031 


1996 Oct 13 


5827 


A 


1.425 


0.171 


0.028 


1996 Oct 14 


5828 


A 


4.860 


0.105 


0.017 



a For previous measurements, see Weiler et al. (1986) and Weiler et al. (1992). 
b The date of the explosion is found to be 1980 October 2 (27 days before optical 
maximum) from the fitting process. This is slightly different from the explosion date 

of 1980 October 1 assumed by Weiler et al. (1986). 

c Measurement presented originally as a < 3a detection which is now interpreted as 

a 3er upper limit. See Footnote |5|. 
d 3cr upper limit. 

e Data from Weiler et al. (1992) that was remeasured for this paper; these values 
supersede those reported in Weiler et al. (1992). See Footnote ||. 
Measurements from Lacey (1997); details in Lacey et al. (1997). 
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Table 3. Fitting Parameters for SN 1980K a 



Parameter 


Value 


Uncertainty Range b 


K x (mJy) 


115 


81-168 


a 


-0.60 


-(0.56-0.67) 


P 


-0.73 


-(0.78-0.68) 


K 2 


1.42 x 10 5 


(0.60-15) xlO 5 


5 


-2.68 


-(2.12-3.11) 


to 


1980 Oct 2 


1980 Sep 11- Nov 1 


M C (M yr- 1 ) 


1.9 x 10~ 5 




X 2 /DOF 


2.57 





a Only data through 1989 July 23 (day 3188) are 
used in the fitting procedure. 

b The error estimates for the parameter values are 

determined using the bootstrap method (see §4). 
c From radio observations prior to 1989 July 23, 

employing Equation 16 from Weiler et al. (1986), 

and assumptions of T = 2.0 x 10 4 K,f shoc k = 1-3 x 

10 4 km s _1 , and w w - md = 10 km s _1 . 



